00010	c	This program computes the matrix and solves for the normalized
00020	c 	age distribution of Sweden.
00030		double precision a(40,30),b(15,15),f(40,15),p(40,15),t(15,15)
00040		double precision x(15),y(15),sgm,s(15,15),z(15),ya(15),za(15)
00050		double precision amax,bmax,ca(15,15),d(15,15),h(15,15)
00060		integer is(15)
00090		n=15
00100		m=15
00110		p(1,1)=140396
00120		p(1,2)=95609
00130		p(1,3)=101699
00140		p(1,4)=96902
00150		p(1,5)=97614
00160		p(1,6)=91277
00170		p(1,7)=81005
00180		p(1,8)=68892
00190		p(1,9)=64301
00200		p(1,10)=59300
00210		p(1,11)=53706
00220		p(1,12)=44752
00230		p(1,13)=38082
00240		p(1,14)=25544
00250		f(1,4)=1653
00260		f(1,5)=10245
00270		f(1,6)=18829
00280		f(1,7)=19237
00290		f(1,8)=14000
00300		f(1,9)=7020
00310		f(1,10)=1621
00320		f(1,15)=72605
00330		p(2,1)=135066
00340		p(2,2)=104580
00350		p(2,3)=97346
00360		p(2,4)=95950
00370		p(2,5)=100288
00380		p(2,6)=90980
00390		p(2,7)=86104
00400		p(2,8)=73275
00410		p(2,9)=66156
00420		p(2,10)=57049
00430		p(2,11)=55371
00440		p(2,12)=46547
00450		p(2,13)=41249
00460		p(2,14)=28211
00470		f(2,4)=1543
00480		f(2,5)=9596
00490		f(2,6)=17431
00500		f(2,7)=18528
00510		f(2,8)=13225
00520		f(2,9)=6305
00530		f(2,10)=1240
00540		f(2,15)=67868
00550		p(3,1)=136048
00560		p(3,2)=103783
00570		p(3,3)=104016
00580		p(3,4)=92228
00590		p(3,5)=95915
00600		p(3,6)=92418
00610		p(3,7)=84428
00620		p(3,8)=77390
00630		p(3,9)=69288
00640		p(3,10)=57836
00650		p(3,11)=52770
00660		p(3,12)=47951
00670		p(3,13)=39872
00680		p(3,14)=29988
00690		f(3,4)=1499
00700		f(3,5)=9737
00710		f(3,6)=18102
00720		f(3,7)=19708
00730		f(3,8)=15073
00740		f(3,9)=7169
00750		f(3,10)=1417
00760		f(3,11)=1
00770		f(3,15)=72706-1
00780		p(4,1)=149560
00790		p(4,2)=110554
00800		p(4,3)=107281
00810		p(4,4)=100565
00820		p(4,5)=96044
00830		p(4,6)=92687
00840		p(4,7)=87918
00850		p(4,8)=79531
00860		p(4,9)=74996
00870		p(4,10)=62790
00880		p(4,11)=54283
00890		p(4,12)=45887
00900		p(4,13)=42391
00910		p(4,14)=31037
00920		f(4,4)=1509
00930		f(4,5)=10043
00940		f(4,6)=19317
00950		f(4,7)=21259
00960		f(4,8)=15652
00970		f(4,9)=7995
00980		f(4,10)=1471
00990		f(4,15)=77246
01000		p(5,1)=149440
01010		p(5,2)=121236
01020		p(5,3)=111919
01030		p(5,4)=103482
01040		p(5,5)=103592
01050		p(5,6)=91662
01060		p(5,7)=88793
01070		p(5,8)=82520
01080		p(5,9)=76806
01090		p(5,10)=69612
01100		p(5,11)=60223
01110		p(5,12)=47745
01120		p(5,13)=40937
01130		p(5,14)=32708
01140		f(5,4)=1465
01150		f(5,5)=10173
01160		f(5,6)=17564
01170		f(5,7)=19757
01180		f(5,8)=15404
01190		f(5,9)=7534
01200		f(5,10)=1506
01210		f(5,15)=73403
01220		p(6,1)=151044
01230		p(6,2)=122818
01240		p(6,3)=120022
01250		p(6,4)=107515
01260		p(6,5)=104536
01270		p(6,6)=97604
01280		p(6,7)=85579
01290		p(6,8)=82043
01300		p(6,9)=78632
01310		p(6,10)=69382
01320		p(6,11)=64376
01330		p(6,12)=51225
01340		p(6,13)=41124
01350		p(6,14)=31866
01360		f(6,4)=1379
01370		f(6,5)=10538
01380		f(6,6)=19392
01390		f(6,7)=19310
01400		f(6,8)=15196
01410		f(6,9)=8209
01420		f(6,10)=1565
01430		f(6,11)=24
01440		f(6,15)=75613-24
01450		p(7,1)=140872
01460		p(7,2)=120437
01470		p(7,3)=118132
01480		p(7,4)=113864
01490		p(7,5)=105820
01500		p(7,6)=97126
01510		p(7,7)=90865
01520		p(7,8)=76985
01530		p(7,9)=75988
01540		p(7,10)=69784
01550		p(7,11)=63256
01560		p(7,12)=55237
01570		p(7,13)=43839
01580		p(7,14)=30498
01590		f(7,4)=1459
01600		f(7,5)=10836
01610		f(7,6)=19838
01620		f(7,7)=20840
01630		f(7,8)=14630
01640		f(7,9)=7527
01650		f(7,10)=1494
01660		f(7,15)=76624
01670		p(8,1)=152319
01680		p(8,2)=116947
01690		p(8,3)=118122
01700		p(8,4)=115204
01710		p(8,5)=113975
01720		p(8,6)=100790
01730		p(8,7)=92831
01740		p(8,8)=85594
01750		p(8,9)=73258
01760		p(8,10)=68766
01770		p(8,11)=64747
01780			p(8,12)=55999
01790		p(8,13)=48132
01800		p(8,14)=33975
01810		f(8,4)=1401
01820		f(8,5)=11468
01830		f(8,6)=20517
01840		f(8,7)=21542
01850		f(8,8)=16909
01860		f(8,9)=7745
01870		f(8,10)=1329
01880	           f(8,11)=2
01890		f(8,15)=80913-2
01900		p(9,1)=167077
01910		p(9,2)=130095
01920		p(9,3)=115945
01930		p(9,4)=116772
01940		p(9,5)=114830
01950		p(9,6)=110321
01960		p(9,7)=97272
01970		p(9,8)=87450
01980		p(9,9)=82058
01990		p(9,10)=67026
02000		p(9,11)=64101
02010		p(9,12)=57578
02020		p(9,13)=48085
02030		p(9,14)=37842
02040		f(9,4)=1319
02050		f(9,5)=11917
02060		f(9,6)=22828
02070		f(9,7)=23484
02080		f(9,8)=17825
02090		f(9,9)=9360
02100		f(9,10)=1491
02110		f(9,11)=12
02120		f(9,15)=88236-12
02130		p(10,1)=193789
02140		p(10,2)=144954
02150		p(10,3)=129598
02160		p(10,4)=118385
02170		p(10,5)=114414
02180		p(10,6)=109835
02190		p(10,7)=106580
02200		p(10,8)=91785
02210		p(10,9)=82742
02220		p(10,10)=75948
02230		p(10,11)=62216
02240		p(10,12)=57307
02250		p(10,13)=50173
02260		p(10,14)=38249
02270		f(10,4)=1395
02280		f(10,5)=12843
02290		f(10,6)=24283
02300		f(10,7)=26277
02310		f(10,8)=19333
02320		f(10,9)=9693
02330		f(10,10)=1653
02340		f(10,11)=6
02350		f(10,15)=95483-6
02360		p(11,1)=192484
02370		p(11,2)=168502
02380		p(11,3)=143815
02390		p(11,4)=131469
02400		p(11,5)=110950
02410		p(11,6)=110187
02420		p(11,7)=106221
02430			p(11,8)=99806
02440		p(11,9)=86794
02450		p(11,10)=76561
02460		p(11,11)=69918
02470		p(11,12)=54714
02480		p(11,13)=49097
02490		p(11,14)=39004
02500		f(11,4)=1189
02510		f(11,5)=11616
02520		f(11,6)=23240
02530		f(11,7)=25580
02540		f(11,8)=20422
02550		f(11,9)=9942
02560		f(11,10)=1527
02570		f(11,11)=1
02580		f(11,15)=93516
02590		p(12,1)=195786
02600		p(12,2)=166813
02610		p(12,3)=164847
02620		p(12,4)=145525
02630		p(12,5)=122754
02640		p(12,6)=106474
02650		p(12,7)=105775
02660		p(12,8)=99926
02670		p(12,9)=94657
02680		p(12,10)=79993
02690		p(12,11)=70016
02700		p(12,12)=61335
02710		p(12,13)=47056
02720		p(12,14)=38616
02730		f(12,4)=1222
02740		f(12,5)=12707
02750		f(12,6)=23497
02760		f(12,7)=26274
02770		f(12,8)=21292
02780		f(12,9)=11387
02790		f(12,10)=1652
02800		f(12,15)=98031
02810		p(13,1)=192553
02820		p(13,2)=171175
02830		p(13,3)=163734
02840		p(13,4)=166204
02850		p(13,5)=137315
02860		p(13,6)=119395
02870		p(13,7)=101892
02880		p(13,8)=99702
02890		p(13,9)=94709
02900		p(13,10)=87520
02910		p(13,11)=73863
02920		p(13,12)=61607
02930		p(13,13)=56001
02940		p(13,14)=36456
02950		f(13,4)=1001
02960		f(13,5)=12288
02970		f(13,6)=24381
02980		f(13,7)=24796
02990		f(13,8)=20476
03000		f(13,9)=10756
03010		f(13,10)=1660
03020		f(13,15)=95358
03030		p(14,1)=207417
03040		p(14,2)=172299
03050		p(14,3)=168329
03060		p(14,4)=165942
03070		p(14,5)=158313
03080		p(14,6)=132397
03090		p(14,7)=115119
03100		p(14,8)=97082
03110		p(14,9)=95245
03120		p(14,10)=88732
03130		p(14,11)=81737
03140		p(14,12)=66141
03150		p(14,13)=56998
03160		p(14,14)=42056
03170		f(14,4)=949
03180		f(14,5)=13100
03190		f(14,6)=26471
03200		f(14,7)=27532
03210		f(14,8)=20463
03220		f(14,9)=11145
03230		f(14,10)=1619
03240		f(14,15)=101279
03250		p(15,1)=216315
03260		p(15,2)=183986
03270		p(15,3)=168029
03280		p(15,4)=169635
03290		p(15,5)=157581
03300		p(15,6)=152632
03310		p(15,7)=128452
03320		p(15,8)=110059
03330		p(15,9)=92878
03340		p(15,10)=89144
03350		p(15,11)=82818
03360		p(15,12)=73395
03370		p(15,13)=59771
03380		p(15,14)=42627
03390		f(15,4)=872
03400		f(15,5)=12515
03410		f(15,6)=29167
03420		f(15,7)=30657
03430		f(15,8)=22918
03440		f(15,9)=11079
03450		f(15,10)=1713
03460		f(15,15)=108921
03470		p(16,1)=230484
03480		p(16,2)=192154
03490		p(16,3)=178858
03500		p(16,4)=170580
03510		p(16,5)=160148
03520		p(16,6)=151127
03530		p(16,7)=147806
03540		p(16,8)=122222
03550		p(16,9)=104463
03560		p(16,10)=86745
03570		p(16,11)=81985
03580		p(16,12)=73542
03590		p(16,13)=63654
03600		p(16,14)=45852
03610		f(16,4)=835
03620		f(16,5)=12483
03630		f(16,6)=28308
03640		f(16,7)=34257
03650		f(16,8)=25928
03660		f(16,9)=12657
03670		f(16,10)=1704
03680		f(16,11)=30
03690		f(16,15)=116172
03700		p(17,1)=251575
03710		p(17,2)=200912
03720		p(17,3)=184250
03730		p(17,4)=175990
03740		p(17,5)=158486
03750		p(17,6)=154845
03760		p(17,7)=146071
03770		p(17,8)=142562
03780		p(17,9)=117148
03790		p(17,10)=98878
03800		p(17,11)=80214
03810		p(17,12)=74051
03820		p(17,13)=64086
03830		p(17,14)=52170
03840		f(17,4)=1001
03850		f(17,5)=13581
03860		f(17,6)=30500
03870		f(17,7)=35549
03880		f(17,8)=31700
03890		f(17,9)=15952
03900		f(17,10)=2123
03910		f(17,15)=130406
03920		p(18,1)=269045
03930		p(18,2)=226408
03940		p(18,3)=196233
03950		p(18,4)=180712
03960		p(18,5)=172262
03970		p(18,6)=154018
03980		p(18,7)=149625
03990		p(18,8)=142009
04000		p(18,9)=136814
04010		p(18,10)=112852
04020		p(18,11)=93399
04030		p(18,12)=73152
04040		p(18,13)=65071
04050		p(18,14)=54317
04060		f(18,4)=1202
04070		f(18,5)=15299
04080			f(18,6)=30548
04090		f(18,7)=35794
04100		f(18,8)=30675
04110		f(18,9)=18205
04120		f(18,10)=2354
04130		f(18,15)=134077
04140		p(19,1)=246626
04150		p(19,2)=241114
04160		p(19,3)=216979
04170		p(19,4)=187815
04180		p(19,5)=167808
04190		p(19,6)=159602
04200		p(19,7)=144150
04210		p(19,8)=140470
04220		p(19,9)=132516
04230		p(19,10)=127302
04240		p(19,11)=104090
04250		p(19,12)=84371
04260		p(19,13)=63471
04270		p(19,14)=53156
04280		f(19,4)=1094
04290		f(19,5)=13723
04300		f(19,6)=28669
04310		f(19,7)=31835
04320		f(19,8)=28347
04330		f(19,9)=15394
04340		f(19,10)=2296
04350		f(19,15)=121358
04360		p(20,1)=266485
04370		p(20,2)=227493
04380		p(20,3)=233571
04390		p(20,4)=211538
04400		p(20,5)=177856
04410		p(20,6)=156805
04420		p(20,7)=150254
04430		p(20,8)=136710
04440		p(20,9)=133114
04450		p(20,10)=125299
04460		p(20,11)=119287
04470		p(20,12)=95660
04480			p(20,13)=75041
04490		p(20,14)=53141
04500		f(20,4)=1751
04510		f(20,5)=17865
04520		f(20,6)=32718
04530		f(20,7)=33529
04540		f(20,8)=29492
04550		f(20,9)=17297
04560		f(20,10)=2381
04570		f(20,11)=10
04580		f(20,15)=135033
04590		p(21,1)=279046
04600		p(21,2)=239647
04610		p(21,3)=219509
04620		p(21,4)=225513
04630		p(21,5)=199236
04640		p(21,6)=167738
04650		p(21,7)=149985
04660		p(21,8)=143954
04670		p(21,9)=130593
04680		p(21,10)=126807
04690		p(21,11)=118370
04700		p(21,12)=110737
04710		p(21,13)=85944
04720		p(21,14)=63991
04730		f(21,4)=2171
04740		f(21,5)=20286
04750		f(21,6)=33564
04760		f(21,7)=33772
04770		f(21,8)=27898
04780		f(21,9)=14885
04790		f(21,10)=2310
04800		f(21,15)=134886
04810		p(22,1)=279795
04820		p(22,2)=249473
04830		p(22,3)=226723
04840		p(22,4)=208137
04850		p(22,5)=203778
04860		p(22,6)=178214
04870		p(22,7)=153914
04880		p(22,8)=140601
04890		p(22,9)=135630
04900		p(22,10)=123144
04910		p(22,11)=118812
04920		p(22,12)=109134
04930		p(22,13)=99418
04940		p(22,14)=73090
04950		f(22,4)=2226
04960		f(22,5)=21487
04970		f(22,6)=35496
04980		f(22,7)=34153
04990		f(22,8)=27298
05000		f(22,9)=14922
05010		f(22,10)=2022
05020		f(22,15)=137604
05030		p(23,1)=289076
05040		p(23,2)=256882
05050		p(23,3)=239978
05060		p(23,4)=212072
05070			p(23,5)=182128
05080		p(23,6)=181970
05090		p(23,7)=165601
05100		p(23,8)=145286
05110		p(23,9)=133350
05120		p(23,10)=128800
05130		p(23,11)=116303
05140		p(23,12)=110489
05150		p(23,13)=98713
05160		p(23,14)=85966
05170		f(23,4)=2255
05180		f(23,5)=19346
05190		f(23,6)=34889
05200		f(23,7)=34567
05210		f(23,8)=26447
05220		f(23,9)=13850
05230		f(23,10)=1938
05240		f(23,15)=133292
05250		p(24,1)=282868
05260		p(24,2)=266543
05270		p(24,3)=247075
05280		p(24,4)=226469
05290		p(24,5)=188173
05300		p(24,6)=164408
05310		p(24,7)=169631
05320		p(24,8)=157474
05330		p(24,9)=138457
05340		p(24,10)=127053
05350		p(24,11)=121796
05360		p(24,12)=108355
05370		p(24,13)=100229
05380		p(24,14)=85258
05390		f(24,4)=2813
05400		f(24,5)=20513
05410		f(24,6)=31702
05420		f(24,7)=34814
05430		f(24,8)=27526
05440		f(24,9)=13846
05450		f(24,10)=1723
05460		f(24,15)=132937
05470		p(25,1)=289567
05480		p(25,2)=268139
05490		p(25,3)=260969
05500		p(25,4)=236210
05510		p(25,5)=209290
05520		p(25,6)=177613
05530			p(25,7)=158435
05540		p(25,8)=163857
05550		p(25,9)=151510
05560		p(25,10)=132816
05570		p(25,11)=121001
05580		p(25,12)=114269
05590		p(25,13)=99228
05600		p(25,14)=87837
05610		f(25,4)=3623
05620		f(25,5)=24443
05630		f(25,6)=34422
05640		f(25,7)=31472
05650		f(25,8)=27484
05660		f(25,9)=13974
05670		f(25,10)=1596
05680		f(25,15)=137014
05690		p(26,1)=297079
05700		p(26,2)=273569
05710		p(26,3)=261573
05720		p(26,4)=246657
05730		p(26,5)=213542
05740		p(26,6)=193700
05750		p(26,7)=168831
05760		p(26,8)=152333
05770		p(26,9)=157566
05780		p(26,10)=145263
05790		p(26,11)=126433
05800		p(26,12)=113488
05810		p(26,13)=104480
05820		p(26,14)=87266
05830		f(26,4)=4107
05840		f(26,5)=25206
05850		f(26,6)=36198
05860		f(26,7)=31778
05870		f(26,8)=23447
05880		f(26,9)=13252
05890		f(26,10)=1539
05900		f(26,15)=135527
05910		p(27,1)=303982
05920		p(27,2)=285787
05930		p(27,3)=269113
05940		p(27,4)=250907
05950		p(27,5)=229124
05960		p(27,6)=202859
05970		p(27,7)=186612
05980		p(27,8)=163048
05990		p(27,9)=146989
06000		p(27,10)=151584
06010		p(27,11)=138660
06020		p(27,12)=119059
06030		p(27,13)=104545
06040		p(27,14)=92267
06050		f(27,4)=4840
06060		f(27,5)=27318
06070		f(27,6)=36056
06080		f(27,7)=32223
06090		f(27,8)=23159
06100		f(27,9)=10946
06110		f(27,10)=1421
06120		f(27,15)=135963
06130		p(28,1)=293099
06140		p(28,2)=293356
06150		p(28,3)=280973
06160		p(28,4)=260955
06170		p(28,5)=236085
06180		p(28,6)=218687
06190		p(28,7)=195927
06200		p(28,8)=180723
06210		p(28,9)=157539
06220		p(28,10)=141704
06230		p(28,11)=145093
06240		p(28,12)=130997
06250		p(28,13)=109873
06260		p(28,14)=92311
06270		f(28,4)=4740
06280		f(28,5)=25218
06290		f(28,6)=33427
06300		f(28,7)=28963
06310		f(28,8)=21479
06320		f(28,9)=10138
06330		f(28,10)=1066
06340		f(28,15)=125031
06350		p(29,1)=272778
06360		p(29,2)=283794
06370		p(29,3)=289388
06380		p(29,4)=273864
06390		p(29,5)=249859
06400		p(29,6)=227216
06410		p(29,7)=211243
06420		p(29,8)=189530
06430		p(29,9)=174611
06440		p(29,10)=151740
06450		p(29,11)=135632
06460		p(29,12)=137058
06470		p(29,13)=120805
06480		p(29,14)=97405
06490		f(29,4)=4892
06500		f(29,5)=26459
06510		f(29,6)=33050
06520		f(29,7)=28264
06530		f(29,8)=20114
06540		f(29,9)=9489
06550		f(29,10)=1043
06560		f(29,15)=123311
06570		p(30,1)=268186
06580		p(30,2)=266013
06590		p(30,3)=280291
06600		p(30,4)=283828
06610		p(30,5)=263443
06620		p(30,6)=240861
06630		p(30,7)=220927
06640		p(30,8)=205949
06650		p(30,9)=184061
06660		p(30,10)=168952
06670		p(30,11)=145697
06680		p(30,12)=128394
06690		p(30,13)=126933
06700		p(30,14)=107459
06710		f(30,4)=4788
06720		f(30,5)=23093
06730		f(30,6)=28573
06740		f(30,7)=23732
06750		f(30,8)=16937
06760		f(30,9)=7760
06770		f(30,10)=869
06780		f(30,15)=105752
06790		p(31,1)=225360
06800		p(31,2)=263570
06810		p(31,3)=263732
06820		p(31,4)=275986
06830		p(31,5)=274636
06840		p(31,6)=254652
06850		p(31,7)=233884
06860		p(31,8)=215158
06870		p(31,9)=200527
06880		p(31,10)=178384
06890		p(31,11)=162261
06900		p(31,12)=138091
06910		p(31,13)=119074
06920		p(31,14)=112880
06930		f(31,4)=4965
06940		f(31,5)=21844
06950		f(31,6)=25340
06960		f(31,7)=20582
06970		f(31,8)=13647
06980		f(31,9)=6124
06990		f(31,10)=654
07000		f(31,15)=93156
07010		p(32,1)=203827
07020		p(32,2)=222951
07030		p(32,3)=263269
07040		p(32,4)=262311
07050		p(32,5)=271887
07060		p(32,6)=271172
07070		p(32,7)=252357
07080		p(32,8)=231055
07090		p(32,9)=211340
07100		p(32,10)=195464
07110		p(32,11)=172335
07120		p(32,12)=154422
07130		p(32,13)=128543
07140		p(32,14)=105923
07150		f(32,4)=4599
07160		f(32,5)=20406
07170		f(32,6)=24712
07180		f(32,7)=19753
07190		f(32,8)=12189
07200		f(32,9)=4902
07210		f(32,10)=502
07220		f(32,15)=87063
07230		p(33,1)=211376
07240		p(33,2)=202251
07250		p(33,3)=222285
07260		p(33,4)=261971
07270		p(33,5)=260442
07280		p(33,6)=269902
07290		p(33,7)=268777
07300		p(33,8)=249416
07310		p(33,9)=227559
07320		p(33,10)=206651
07330		p(33,11)=189117
07340		p(33,12)=164479
07350		p(33,13)=143835
07360		p(33,14)=114510
07370		f(33,4)=5436
07380		f(33,5)=23456
07390		f(33,6)=29012
07400		f(33,7)=23267
07410		f(33,8)=13718
07420		f(33,9)=4837
07430		f(33,10)=430
07440		f(33,15)=100156
07450		p(34,1)=277700
07460		p(34,2)=217359
07470		p(34,3)=203335
07480		p(34,4)=222289
07490		p(34,5)=260561
07500		p(34,6)=260021
07510		p(34,7)=269972
07520		p(34,8)=267645
07530		p(34,9)=247440
07540		p(34,10)=224942
07550		p(34,11)=200157
07560		p(34,12)=180259
07570		p(34,13)=153415
07580		p(34,14)=129587
07590		f(34,4)=6886
07600		f(34,5)=31650
07610		f(34,6)=37682
07620		f(34,7)=30344
07630		f(34,8)=18423
07640		f(34,9)=6002
07650		f(34,10)=438
07660		f(34,15)=131425
07670		p(35,1)=301458
07680		p(35,2)=278985
07690		p(35,3)=218506
07700		p(35,4)=206242
07710		p(35,5)=230626
07720		p(35,6)=269656
07730		p(35,7)=263285
07740		p(35,8)=271395
07750		p(35,9)=267140
07760		p(35,10)=245173
07770		p(35,11)=219875
07780		p(35,12)=195075
07790		p(35,13)=172055
07800		p(35,14)=140628
07810		f(35,4)=7814
07820		f(35,5)=29060
07830		f(35,6)=34695
07840		f(35,7)=24946
07850		f(35,8)=14768
07860		f(35,9)=5071
07870		f(35,10)=391
07880		f(35,15)=116745
07890		p(36,1)=260125
07900		p(36,2)=301271
07910		p(36,3)=278977
07920		p(36,4)=220918
07930		p(36,5)=215298
07940		p(36,6)=235544
07950		p(36,7)=271019
07960		p(36,8)=262533
07970		p(36,9)=268811
07980		p(36,10)=264178
07990		p(36,11)=240347
08000		p(36,12)=212878
08010		p(36,13)=185032
08020		p(36,14)=157269
08030		f(36,4)=8419
08040		f(36,5)=28513
08050		f(36,6)=31197
08060		f(36,7)=23030
08070		f(36,8)=12159
08080		f(36,9)=3916
08090		f(36,10)=298
08100		f(36,15)=107532
08110		p(37,1)=253236
08120		p(37,2)=260407
08130		p(37,3)=301498
08140		p(37,4)=281815
08150		p(37,5)=227717
08160		p(37,6)=217682
08170		p(37,7)=236238
08180		p(37,8)=270408
08190		p(37,9)=260812
08200		p(37,10)=265861
08210		p(37,11)=259584
08220		p(37,12)=233660
08230		p(37,13)=203042
08240		p(37,14)=170792
08250		f(37,4)=10288
08260		f(37,5)=30125
08270		f(37,6)=30419
08280		f(37,7)=19980
08290		f(37,8)=10672
08300		f(37,9)=3138
08310		f(37,10)=227
08320		f(37,15)=104849
08330		p(38,1)=270162
08340		p(38,2)=255084
08350		p(38,3)=262124
08360		p(38,4)=306081
08370		p(38,5)=288935
08380		p(38,6)=230884
08390		p(38,7)=218454
08400		p(38,8)=236060
08410		p(38,9)=269534
08420		p(38,10)=259637
08430		p(38,11)=261750
08440		p(38,12)=252959
08450		p(38,13)=224007
08460		p(38,14)=188863
08470		f(38,4)=14881
08480		f(38,5)=40711
08490		f(38,6)=35552
08500		f(38,7)=19515
08510		f(38,8)=9275
08520		f(38,9)=2666
08530		f(38,10)=185
08540		f(38,15)=122785
08550		nn=36
08560		n1=1         
08570		sgm=0
08580		do 3000 j=1,10
08590	3000	sgm=sgm+p(nn+1,j)
08600		do 2000 i=1,10
08610	2000	x(i)=p(nn+1,i)/sgm
08620		do 100 i=n1,nn
08630		do 200 j=1,13
08640	200	a(i,14+j)=p(i+1,j+1)/p(i,j)
08650		do 300 j=1,14
08660	300	a(i,j)=f(i,j)*p(i+1,1)/(f(i,15)*p(i,j))
08670		goto 100
08680		write(3,1250)(a(i,j),j=1,14)
08690		write(3,1250)(a(i,j),j=15,28)
08700		write(3,1280)
08710		write(3,1280)
08720		write(3,1280)
08730		write(3,1280)
08740	100	continue
08750		do 400 i=1,14
08760		s(1,i)=a(nn,i)
08770		if(i .eq.14) goto 400
08780		s(i+1,i)=a(nn,i+14)
08790	400	continue
08800		n1=nn-n1  
08810		na=nn+1
08820		do 700 i=1,n1
08822		do 9 j=1,15
08824	9	is(j)=j
08830		n2=nn-i+1
08840		write(3,1280)
08850		do 80 j=1,10
08860		do 90 k=1,10
08870	90	t(k,j)=s(k,j)
08880	80	continue
08890		do 800 j=1,9
08900		do 900 k=1,10
08910	900	s(k,j)=t(k,1)*a(nn-i,j)+t(k,j+1)*a(nn-i,14+j)
08920	800	continue
08930		do 1000 k=1,10
08940		s(k,10)=a(nn-i,10)*t(k,1)
08950	1000	continue
08960	8	n5=n2-1
08970		write(3,1290)na,n5
08980		write(3,1280)
08990	1290	format(1h ,60x,'h(',i2,',',i2,')')
09000	450	do 500 k=1,10
09010		do 600 j=1,10
09020	600	t(k,j)=s(k,j)
09030	500	write(3,1260)(s(k,j),j=1,10)
09040		write(3,1280)
09050		icn=0
09060	4005	lc=1
09070		do 4010 k1=1,10
09080	4010	za(k1)=.1
09090	4800	sgm=0
09100		do 4050 l=1,10
09110		y(l)=0
09120		do 4100 j=1,10
09130	4100	y(l)=y(l)+t(l,j)*za(j)
09140	4050	sgm=sgm+y(l)
09150		do 4500 j=1,10
09160		y(j)=y(j)/sgm
09170	4500	z(j)=(x(j)-y(j))/x(j)
09180		sgm=0
09190		do 4600 j=1,10
09200	4600	sgm=sgm+abs(za(j)-y(j))
09210		if(sgm .le. .00001) goto 4450
09220		lc=lc+1
09230		if(lc .gt. 30) goto 4450
09240		do 4700 j=1,10
09250	4700	za(j)=y(j)
09260		goto 4800
09270	4450	do 4400 j=1,10
09280	4400	write(3,1270)x(j),y(j),z(j)
09290	1285	format(1h ,5(5x,f15.8))
09300		write(3,1300)lc
09310	1300	format(1h ,10x,i3)
09320		write(3,1280)
09330		icn=icn+1
09340		if(icn .eq. 2) goto 715
09350		do 710 i1=1,10
09360		do 720 j1=1,10
09370		t(i1,j1)=0
09380		do 730 k1=1,10
09390	730	t(i1,j1)=t(i1,j1)+s(i1,k1)*s(j1,k1)
09400	720	continue
09410	710	continue
09420		goto 4005
09430	715	sgm=0
09440		do 740 i1=1,10
09450		z(i1)=0
09460		do 750 j1=1,10
09470	750	z(i1)=z(i1)+s(j1,i1)
09480		do 760 j1=1,10
09490	760	b(j1,i1)=s(j1,i1)/z(i1)
09500	740	continue
09510		write(3,1280)
09520		do 5 j=1,9
09530		do 10 k=1,10
09540	10	t(j+1,k)=b(k,1)-b(k,j+1)
09550		t(j+1,11)=0
09560		do 15 k=1,10
09570	15	t(j+1,11)=t(j+1,11)+(b(k,1)-b(k,j+1))*(b(k,1)+b(k,j+1))/2
09580		write(3,1261)(t(j+1,k),k=1,11)
09590	1261	format(1h ,10(1x,f11.8),1x,f11.6)
09610		za(j+1)=t(j+1,11)
09615		do 999 k=1,11
09620	999	h(j+1,k)=t(j+1,k)
09630	5	continue
09640		za(1)=1
09650		write(3,1280)
09660		do 20 j=1,11
09670		h(1,j)=1
09680	20	t(1,j)=1
09682		do 21 j=1,10
09684		do 22 k=1,10
09686	22	ca(j,k)=0
09688	21	continue
09690		np=10
09700		m1=11
09710		ierr=0
09720		call detgau(np,m1,t,ierr,amax,bmax,ca,is)
09730		if(ierr .eq. 1) goto 700
09735		call bksub(np,m1,t)
09738		write(3,1280)
09740		write(3,1260)(t(k,11),k=1,10)
09750		ica=0
09760		do 31 k=1,10
09770	31	y(k)=t(k,11)
09780	14	do 3 j=1,10
09785		sgm=0
09790		do 4 l=1,10
09800	4	sgm=sgm+h(j,l)*y(l)
09820	3	ya(j)=za(j)-sgm
09830		do 11 j=1,10
09840	11	ca(j,11)=ya(is(j))
09850		call fwsub(np,m1,ca)
09860		do 12 j=1,10
09870	12	t(j,11)=ca(j,11)
09880		call bksub(np,m1,t)
09890		sgm=0
09900		do 13 j=1,10
09910		sgm=sgm+abs(t(j,11))
09920	13	y(j)=y(j)+t(j,11)
09930		ica=ica+1
09940		if(ica .gt. 30) goto 18
09950		if(sgm .lt. .000001) goto 18
09952		sgm=0
09953		do 26 j=1,10
09954	26	sgm=sgm+y(j)
09955		do 27 j=1,10
09956	27	y(j)=y(j)/sgm
09957		goto 14
09960	18	sgm=0
09970		do 16 j=1,10
09980	16	sgm=sgm+y(j)
09990		do 17 j=1,10
10000	17	y(j)=y(j)/sgm
10010		do 30 j=1,10
10020	30	z(j)=(x(j)-y(j))/x(j)
10030		do 35 j=1,10
10040	35	write(3,1270)x(j),y(j),z(j)
10042		write(3,1270)sgm,sgm,sgm
10045		write(3,1300) ica
10050		write(3,1280)
10060		write(3,1280)
10070		write(3,1280)
10080		write(3,1280)
10090	700	continue
10100		goto 8000
10110		n3=1
10120		write(3,1290)na,n3
10130		do 1200 i=1,10
10140		write(3,1260)(s(i,j),j=1,10)
10150	1200	continue
10160		write(3,1280)
10170	1250	format(1h ,14(1x,f8.5))
10180	1260	format(1h ,10(1x,e12.6))
10190	1270	format(1h ,3(5x,e15.6))
10200	1280	format(1h ,5x)
10210		write(3,1280)
10220		do 4000 i=1,10
10230	4000	z(i)=.1
10240	    	do 4550 i=1,10
10250		y(i)=0
10260		do 5000 j=1,10
10270	5000	y(i)=y(i)+s(i,j)*z(j)
10280	4550	continue
10290		sgm=0
10300		do 5500 i=1,10
10310	5500	sgm=sgm+y(i)
10320		do 6000 i=1,10
10330		y(i)=y(i)/sgm
10340	6000	write(3,1270) y(i)
10350		write(3,1280)
10360		sgm=0
10370		do 6100 i=1,10
10380		ya(i)=s(i,5)
10390	6100	sgm=sgm+ya(i)
10400		do 6200 i=1,10
10410		z(i)=(x(i)-y(i))/x(i)
10420		ya(i)=ya(i)/sgm
10430		za(i)=(x(i)-ya(i))/x(i)
10440	6200	write(3,1285)x(i),y(i),z(i),ya(i),za(i)
10450	8000	stop
10460		end
10470	
10480		subroutine detgau(n,m,t,ierr,amax,bmax,ca,is)
10490		double precision amax,f,c,bmax
10500		double precision t(15,15),ca(15,15)
10510		integer is(15)
10520		det=1
10530		na=n-1
10540		do 4 k=1,na
10550		amax=abs(t(k,k))
10560		ic=k
10570		do 3 i=k,n
10580		if(abs(t(k,i)) .le. amax) goto 3
10590		amax=abs(t(k,i))
10600	3	continue
10610		amax=abs(t(k,k))/amax
10620		kp1=k+1
10630		do 5 i=kp1,n
10640		bmax=abs(t(i,k))
10650		do 6 j=kp1,n
10660		if(abs(t(i,j)) .le. bmax) goto 6
10670		bmax=abs(t(i,j))
10680	6	continue
10690		if(abs(t(i,k)/bmax) .le. amax) goto 5
10700		amax=abs(t(i,k)/bmax)
10710		ic=i
10720	5	continue
10730		if(ic .eq. k) goto 9
10740		l=is(ic)
10750		is(ic)=is(k)
10760		is(k)=l
10770	7	do 8 j=1,m
10772		f=ca(k,j)
10774		ca(k,j)=ca(ic,j)
10776		ca(ic,j)=f
10780		c=t(k,j)
10790		t(k,j)=t(ic,j)
10800	8	t(ic,j)=c
10810	9	amax=t(k,k)
10830		do 1 i=kp1,n
10840	1	ca(i,k)=t(i,k)/t(k,k)
10850	 	ca(k,k)=1 
10855		if(abs(amax) .lt. .000000001) goto 4
10860		do 11 i=kp1,n
10870		if(t(i,k) .eq. 0) goto 11
10880		do 12 j=kp1,m
10890	12	t(i,j)= t(i,j)-t(k,j)*t(i,k)/amax
10900		t(i,k)=0
10910	11	continue
10920	4	continue
10925		ca(10,10)=1
10930		return
10940		end
10950		
10960		subroutine bksub(n,m,t)
10970		double precision t(15,15),fac
10980	   	np1=n+1
10990		do 15 k=np1,m
11000		do 17 i=1,n
11010		ii=n-i+1
11020		fac=0
11030		if(ii .eq. n) goto 17
11040		iip1=ii+1
11050		do 16 j=iip1,n
11060	16	fac=fac+t(ii,j)*t(j,k)
11070	17	t(ii,k)=(t(ii,k)-fac)/t(ii,ii)
11080	15	continue
11090	20	return
11100		end
11110	
11120		subroutine fwsub(n,m,t)
11130		double precision t(15,15),fac
11140		np1=n+1
11150		do 10 k=np1,m
11160		do 20 i=1,n
11170		fac=0
11180		if(i .eq. 1) goto 20
11190		ii=i-1
11200		do 30 j=1,ii
11210	30	fac=fac+t(i,j)*t(j,k)
11220	20	t(i,k)=(t(i,k)-fac)/t(i,i)
11230	10	continue
11240		return
11250		end